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Abstract 

The binding energies of di- hadronic states have been calculated assuming a 'molecular' 
interaction provided by the asymptotic expression of the residual confined gluon exchange 
potential between the component hadrons in the system. Meson- meson and meson- baryon 
states have been studied in detail and a mass formula has been used to calculate total mass 
of the 'molecules'. 
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I. Introduction 

Quantum chromodynamics allows for the existence of various multi- quark states, and the 
same can also be predicted by various phenomenological models [1-6]. While the properties 
of mesons(qg) and baryons(qqq) have been well documented both in theory and experiment, 
the exotics ( states that do not fit into qq or qqq states ) remain lesser known and understood. 
The exotics include the tetra- quark, penta- quark and hexa- quark states, hybrids such as 
qqg and gg or ggg glueball states [7]. The recent reports on discovery of the pentaquark 
@ + (1540) [8-10] have revived interest in the study of di- hadronic states in particle physics. 
The baryon + which goes to K + n cannot be a 3-quark state and its minimal quark content 
is uudds, making it manifestly exotic. 

Investigations into the existence of multi- quark states using the quark and bag models 
[11,12] began in the early days of QCD. Non- relativistic potential models have also been 
used to study these systems [2,13]. Recently tetraquark states have been studied extensively, 
specially with regard to the X(3872) particle [14-16]. QCD spectral sum rules have been used 
to test the nature of the heavy meson. The tetraquark has been studied in the framework 
of the string(flux tube) model [17] where it is assumed that di- quark- di-antiquark are 
connected by colour flux. Heavy chiral unitary approach has been used to investigate inter- 
actions between heavy vector mesons and light pseudoscalar mesons [18]. S-wave tetraquarks 
have also been studied using chromomagnetic interaction with full account for flavour sym- 
metry breaking [19]. The spectra of tetraquark and pentaquark states has been evaluated. 
Complete classification has been provided in terms of spin flavour, color and spatial degrees 
of freedom for tetraquarks [20] The same has been done for 5 quark systems in terms of spin 
flavour SU(6) representation [21]. However, little success has been achieved in understanding 
tetra- quark and penta- quark states as di- hadronic molecules due to the non-perturbative 
nature of QCD at the hadronic scale. 

The present paper seeks to examine the tetraquark and pentaquark states as di- hadronic 
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molecules (meson- meson and meson- baryon). A molecular interaction of the Van der Waal's 
type is assumed between the constituent hadrons [22,23]. Binding energy of the molecules is 
calculated using two different wavefunctions obtained from the Statistical model [24]. Spin 
hyperfine interactions have also been taken into consideration, and a mass formula employed 
to calculate total mass of the molecule. Results obtained by using the two wavefunctions are 
compared to each other and to available experimental data. The approach has been clearly 
outlined in 'Method', with the calculations being presented in the 'Results' section. The 
work has been rounded off in 'Conclusion' which comprises a detailed analysis of the nature 
of results obtained and how well they compare to experimental data and results obtained by 
other groups. 



II. Method 



The di- hadron is constituted of either a meson- meson or a meson- baryon held together 
by a Van der Waal type of interaction. The formula for calculating the low- lying di- hadronic 
molecular mass is taken as 

M Tot = M 1 + M 2 + E BE + E SD (1) 

where Mi and M 2 represent the masses of the constituent hadrons respectively, and Ebe rep- 
resents the binding energy of the di- hadronic system and Esd represents the spin- dependent 
term, taken separately. 

Ebe can be calculated using the formula 

E BE =<y{r 12 )\V{r l2 )\y{r l2 )> (2) 

where r 12 is taken to be the radius parameter of the di- hadronic molecule and V(r 12 ) is 
the di- hadronic molecular potential [22, 23]. ^(r^) is the wavefunction of the di- hadronic 
state [24]. The wavefunctions for the ground state of the molecule are obtained from the 
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statistical model as 



and 



l^(n 2 )| 2 = -^{r, - r 12 )* 2 0(r o - r 12 ) (3) 
o47rr ' 



I*(ri2)| 2 = -je(r 2 -r? 2 ) 3/2 ^(ro-r 12 ) (4) 

where r is the radius of the molecule and 6(r — r 12 ) represents the step function. The 
residual interaction of the confined gluon is considered similar to the Van der Waal's interac- 
tion and is assumed as that due to asymptotic expression(ri 2 — > oo) of the residual confined 
one gluon exchange interaction with strength k mo i [22,23]. The potential is given as 

V{ r 12 ) = zh^ e -c*r> 2 /2 (5) 

where k mo i is the residual strength of the strong interaction molecular coupling and C is the 
effective colour screening of the confined gluons. The radius r of the di- hadronic molecule in 
Eqs. 3 and 4 is obtained by employing an additive rule to the radii of constituent hadrons. We 
take r Q = r\ + r 2 , where ri and r 2 represent the individual radii of the hadrons constituting 
the molecule, respectively. 

Now, using Eqs. 2, 3 and 5 we get an integral for E BE which yields 

31 ^>k 

Ebe = —ttt — 0.114286#hpPFQ[(3/2, 1), (2.75, 2.25), -(3\ifRe(5 > (6) 
16r 

where Hyp PFQ is a hypergeometric function. Now, using Eqs. 2,4 and 5 we get 

32*^-0.75 + 0.5/3 e'? O.QQ4Q7e^Erftyp} ^ 
Ebe-— I j 2 w j (7) 

where Erfi[] represents an error function and (5— — ^ and C=50 MeV [25] has been sub- 
stituted. The experimental mass for + is 1540 MeV [8] being known, the values of masses of 
the constituent hadrons (K + and n) are substituted in Eq.l to find the corresponding binding 
energy. This is then used in Eq.7 to calculate the value of k mo i. It is found k moi =0.47, which 
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is fixed as the strength of the residual confined gluonic interaction corresponding to the di- 
hadronic wavefunction in Eq.3. Working the same way, k mo i=0.52 is calculated for the wave- 
function in Eq. 4. The radii considered are r(vr) = 5.38Ge\/ _1 , r(K) = 4.77Ge\/~ 1 , r(p) = 
A.lhGeV- 1 [26], r(D s ) = 4.3Ge\/- 1 [27], r(co) = 4.7Q5GeV-\ r(0) = S.OOGeV- 1 [28] .Sim- 
ilarly r{p) = 6GeV-\ R{n) = 4.7Ge\/- 1 [29], r(A+) = 5.727GeV-\ r(E+) = 3.38QGeV-\ 
r(S°) = 2.404Ge\/- 1 [30], r(E) = S^GeV' 1 [31]. Experimental values of the respective 
meson and baryon masses have been considered for calculation [32]. 

Esd, the contribution to total mass from spin hyperfine interaction is given by the fol- 
lowing expression [33] 



where mi and m,2 are the individual masses of the constituent hadrons in the di hadronic 
molecule, a s is the strong interaction constant and Si and S2 are the spins of the hadrons 
involved. [^(O)] 2 is the value of the di- hadronic wavefunction at the origin. The spin 
contribution on calculation is found to be extremely small in comparison with binding energy 
of the molecules. The values of binding energy and spin- dependent term and total mass 
calculated for the di- hadronic molecules are displayed in Tablel (Tetra- quark states) and 
Table II(Penta- quark states)and Table III (hexa-quark states). Available experimental data 
are presented alongside in the tables to provide a ready comparison with obtained values. 

IV. Conclusion 

The binding energies and masses of exotic states (tetraquark and pentaquark) have been 
calculated by considering them as di-hadronic molecules. The values are calculated using 
two separate wavefunctions (Eq. 3 and Eq. 4) obtained using linear and harmonic potentials 
respectively. The values obtained in both cases are found to be in close agreement with no 
serious discrepancy. This may be considered to suggest that for any wavefunction arising 
from a smooth background potential, the calculated values will be independent of the actual 
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potential. The experimental states are exotics whose spin- parity do not match with the 
expected quark- antiquark structure for mesons or the 3- quark structure for baryons, Hence 
pseudoscalar- pseudoscalar di mesonic states and vector- vector di mesonic states are assigned 
parity charge- conjugation PC '++' and pseudoscalar- vector di mesonic states have PC '+- 
'. The values of J PC have also been indicated in the tables and those experimental states 
whose J PC values match the predicted values have been considered for comparison. In case 
of di-mesonic states, a number of experimentally observed states were available, to which 
the calculated values have been compared, wherever possible. The calculated results are 
found to be in good agreement with experiment. The n — K state is found to have a mass of 
approximately 0.73 GeV which is nearly the same as that calculated by Rai et al(~ 0.72GeV) 
[25]. The binding energies of tetra- quark states are found to lie within a range of 0.094 GeV 
to 0.114 GeV, whilst the masses range from 0.373 to 2.575 GeV. The low- lying pentaquark 
states have been treated as meson- baryon molecules. The experimentally observed + 
state [8] has been taken as reference for calculating the value of the interaction strength 
parameter. The obtained values have been compared to the few available experimental data 
and are found to be in favourable agreement. The calculated total mass for the S — K state 
(1.77 GeV) is found to be close to the estimate given by Rai et al [25]. The well known 
A particle (1.405 GeV) [34] is found to almost exactly match the predicted n — E molecule 
(1.41GeV). The binding energies of meson- baryon molecules are found to range from 0.0848 
to 0.145 GeV whilst the total mass ranges from 1.41 to 3.3 GeV. In calculating the binding 
energies and total mass of hexa- quark particles experimental results were not available for 
comparison. Despite theory suggesting that the hexa- quark state should be more stable than 
the penta- quark state, experimental evidence has been found in support of the latter, not 
the former. Binding energies of the hexa- quark states are found to range from 0.0938GeV 
to 0.133GeV, and mass varies from 2.012GeV to 2.479GeV. 

We have recently successfully predicted the mass of the pentaquark 6 + using a diquark 
approach[35]. This suggests that the same approach may be extended to studying other 
multiquark states as has been done here and we look forward to doing the same in our 
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future works. The good agreement of the theoretically calculated results with available 
experimental data and other theoretical works encourages us to suggest that a number of 
the predicted states may be experimentally detected in the near future. 
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V. Results 



Tablel: Binding Energies and Masses of Di-Meson States 







Linear 


Potential 


Harmonic 


Potential 




States 


jPC 


Ebe + Esd 


M 


Ebe + Esd 


M 


MExp 






(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


n — n 


0++ 


0.0948 


0.3738 


0.0944 


0.3734 


0.4-1.2 


n-K 


0++ 


0.1010 


0.7335 


0.1005 


0.7355 




71 — UJ 


1+- 


0.101 


1.023 


0.1005 


1.0227 




7T — 


1+- 


0.0985 


1.257 


0.0978 


1.256 


1.235 (6i) 


n-D s 


0++ 


0.1059 


2.213 


0.1055 


2.212 


2.317 od: 


71 — p 


1+- 


0.0985 


1.257 


0.0979 


1.256 


1.17 (hj 


K-u 


1+- 


0.1075 


1.3839 


0.1074 


1.3838 


1.4 (tfi) 


K-<j> 


1+- 


0.1049 


1.618 


0.1045 


1.6167 




K-D s 


0++ 


0.1142 


2.575 


0.1132 


2.575 




K-K 


0++ 


0.1075 


1.0948 


0.1074 


1.0947 


1.27 (K x ) 


K- p 


1+- 


0.1076 


1.376 


0.1076 


1.376 


1.41 (K 1} 


P~ P 


1++ 


0.108 


1.659 


0.1078 


1.658 


1.69 (p 3 ) 


p — oj 


2++ 


0.1077 


1.665 


0.1077 


1.665 




P ~ 


2++ 


0.1052 


1.9 


0.1047 


1.898 




0-0 


2++ 


0.1024 


2.1144 


0.102 


2.14 


2-01 (f 2 ) 
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Tablell: Binding Energies and Masses of Di-Baryon States 
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1 01 7 


1 1807 


o 1 m a 


1 1804 
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IT 


_ p 
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1/Z 
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1 1fi79 

J. . J.U I Z 


U .UOO-7 






IT 


_ V 


1/2+ 

1/ z 


D 0855 


1.4118 


0848 

\J • \J UtU 


1 4168 


1 4(Sm) 


/I 




1 /9+ 




9 71 9^ 


1 1 


9 71 9^ 

Z. 1 ±ZU 




/I 


- A + 


l/Z 


n DQ1 8 


9 ^1 04 


n qi if\ 


9 M 70 


9 ^9^ 


7T - 


_ -0 
'- 'c 


3/9+ 




9 741 


n 1 ^9 


9 749 




K 


-p 


1 /9~ 

l/Z 


0Q47 


1 ^9^8 


n oq^ 

U. \JcfO 


1 ^10 




K 


- E 


1/2 - 


0.0904 


1.77 


0899 


1.77 


1.75(E) 


K- 




1/2- 


0.0973 


2.877 


0.0973 


2.877 




K - 


_ ^0 

1 — 'c 


1/2- 


0.145 


3.109 


0.144 


3.054 




P 


-iV 


l/2+,3/2+ 


0.1088 


1.826 


0.1086 


1.823 


1-9(P 35) ) 


P 


-p 


l/2+,3/2+ 


0.0948 


1.8078 


0.0945 


1.8079 




P 


- E 


1/2+ 


0.0898 


2.0508 


0.0902 


2.055 




P~ 


■s+ 


1/2+ 


0.1289 


3.3526 


0.1289 


3.352 




U) 


-TV 


l/2+,3/2+ 


0.1088 


1.83 


0.108 


1.83 




UJ 


- E 


1/2+ 


0.0898 


2.06 


0.089 


2.06 




UJ 


-P 


1/2+3/2+ 


0948 


1.814 


0.095 


1.814 




<f> 


- TV 


l/2+,3/2+ 


0.1055 


2.064 


0.1058 


2.064 







-P 


l/2+,3/2+ 


0.927 


2.05 


0.0928 


2.05 







- E 


1/2+ 


0.088 


2.301 


0.089 


2.301 





Tablelll: Binding Energies and Masses of Di-Baryon States 







Linear 


Potential 


Harmonic 


Potential 




States 


jPC 


Ebe + E S d 


M 


Ebe + E S d 


M 


M Exp 






(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


N-N 




0.133 


2.012 


0.1297 


2.008 




N- E 




0.1106 


2.428 


0.109 


2.241 




P- E 




0.0985 


2.229 


0.0976 


2.229 




E - E 




0.0938 


2.479 


0.0939 


2.479 
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